Matter diffraction at oblique incidence: Higher resolution and the 4 He 3 Efimov state 
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We study the diffraction of atoms and weakly-bound three-atomic molecules from a transmission grating at 
non-normal incidence. Due to the thickness of the grating bars the slits are partially shadowed. Therefore, the 
projected slit width decreases more strongly with the angle of incidence than the projected period, increasing, 
in principle, the experimental resolution. The shadowing, however, requires a revision of the theory of atom 
diffraction. We derive an expression in the style of the Kirchhoff integral of optics and show that the diffraction 
pattern exhibits a characteristic asymmetry which must be accounted for when comparing with experimental 
data. We then analyze the diffraction of weakly bound trimers and show that their finite size manifests itself in 
a further reduction of the slit width by (3/4) (r) where {r} is the average bond length. The improved resolution 
at non-normal incidence may in particular allow to discern, by means of their bond lengths, between the small 
ground state of the helium trimer {(r) « 1 nm, Barletta and Kievsky, Phys. Rev. A 64, 042514 (2001)) and its 
predicted Efimov-type excited state ({r) « 8 nm, ibid.), and in this way to experimentally prove the existence of 
this long-sought Efimov state. 

PACS numbers: 36.40.-c, 03.75.Be, 36.90.+f 



I. INTRODUCTION 

The combination of two unique features makes matter-wave 
diffraction of noble gas trimers an outstanding enterprise. 
Firstly, diffraction presently is the only experimental tech- 
nique which allows to detect such very weakly bound clusters 
and to determine their properties 12> HI • Secondly, the he- 
lium trimer 4 He3 is the only molecule predicted to possess an 
Efimov-type bound state |4|,|5|,|6|] under normal conditions |7]. 

Only recently did diffraction of atomic and molecular 
beams evolve towards a precise experimental technique. Early 
pioneering experiments had proved diffraction for a sodium 
beam through a grating initially fabricated for x-rays |8|], as 
well as for meta-stable helium |9] and for neon 1 10] through 
micrometer double-slits. Later, transmission gratings with a 
period of d — 200 nm had brought an improvement, but fi- 
nally the production of reliable nano-scale transmission grat- 
ings with a period of only d = lOOnm 1 11] paved the way 
for quantitative matter diffraction experiments with an unpar- 
alleled spatial coherence across up to a hundred slits for a 
helium atom beam fl2Tl . Unprecedented, matter diffraction 
also allows to investigate the coherence properties of very 
heavy molecules with many internal degrees of freedom such 
as fullerenes 1 13U14II . 

Classical wave optics can merely serve as an approximation 
to the underlying physical scattering process of atom diffrac- 
tion. The hierarchy of the diffraction peak intensities in a 
diffraction experiment with neutral atoms or molecules is sig- 
nificantly affected by the weak van der Waals surface force, 
which acts on atoms in the vicinity of the material grating. 
This was included in a quantitative theory in Refs. 1151 Hal 
which allowed, in comparison with experimental data, to char- 
acterize this surface force for the noble gases helium, neon, ar- 
gon, and krypton, and the covalently bound D2 molecule 1 12], 
as well as for meta-stable helium and neon I 1 1711 . 

Moreover, there is no analog in wave optics for the diffrac- 
tion of weakly bound small noble gas van der Waals clusters 
such as the helium dimer 4 He2 and trimer 4 He3. Experimen- 



tal evidence for these delicate molecules was for the first time 
unambiguously provided through the mass-selective property 
of grating diffraction |2]. Moreover, the comparatively large 
interatomic distance of (r) = 5.2 nm in 4 He2, implied by the 
small binding energy £\, = -1-1 mK |3], was shown to man- 
ifest itself as an apparent narrowing of the grating slits by 
= (r). It was this size effect which rendered possible the de- 
termination of (r) from experimental diffraction data |3[. This 
allowed for comparison with modern realistic helium-helium 
potentials of Refs. fl^fl9ll . 

Numerical studies of the helium trimer r elying o n these 
helium-helium potentials B H E3 EH E2lS El 
l27l 12811 . quantum Monte Carlo simulations 1 29], as well as 
quantum chemical ab initio calculations l30Tl have long pre- 
dicted two bound states for 4 He3: a ground state at -126 mK 
and a shallow excited state at -2.3 mK |4]. Moreover, the ex- 
cited state is believed to be an Efimov-ty pe s tate. Originally 
in the context of nuclear physics, Efimov 1 3 1 ] had shown that 
if the scattering length of a pair potential exceeds the effec- 
tive range of the potential by far then a universal series of 
bound states exists in the three-body system near the dissoci- 
ation threshold. Examples for such Efimov states, however, 
have been searched for in vain in three-nucleon systems, leav- 
ing the three-atomic helium molecule presently as the only 
candidate. 

The 4 He3 excited state cannot experimentally be distin- 
guished from the ground state by its mass. However, due to 
the large difference in binding energy, both predicted states 
have markedly different interatomic distances: (r) = 0.96 nm 
in the ground state and (r) = 7.97 nm in the excited state |4]. 
Therefore, the size effect, which had previously played the es- 
sential role in the dimer diffraction experiment and which we 
show to be | (r) for a trimer, is expected to render the two 
states distinguishable. Bruch et al have analyzed the mole 
fraction of small helium clusters (including atoms) in a noz- 
zle beam diffraction setup and showed that up to 7% can be 
trimers [32]. This should be an ample amount for a quantita- 
tive analysis. The population ratio of ground state vs. excited 
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state trimers in the beam is, however, not known. It is there- 
fore essential to provide sufficient experimental resolution for 
the small ground state in order to evaluate diffraction data 
from a mixed beam. A limitation in the resolution is posed 
by the period of the grating, d = 100 nm, and the slit width, 
typically so = 60 nm, which are both large compared to the 
ground state size. Transmission gratings with smaller periods 
and slit widths are, however, presently not available. 

To address this issue we consider diffraction from a cus- 
tom transmission grating at oblique (non-normal) incidence, 
i.e. the grating is rotated by an angle 6' about an axis parallel 
to its bars. The grating typically consists of a t as 120 nm thick 
layer of silicon nitride into which the slits have been etched in 
a lithographic production process full . As the grating is ro- 
tated the upstream edges of its bars cut into the incident beam, 
partially shadowing the slits. By this effect the projected 
slit s± width (perpendicular to the incident beam) decreases 
more strongly with 8' than the projected period d ± . Since 
the diffraction pattern is roughly governed by the squared slit 
function 13311 

r i2 

sm(nns ± /d ± ) 
(nns_ L /d±) 

where n denotes the diffraction order, one expects that at non- 
normal incidence more atoms or molecules are diffracted into 
higher diffraction orders. For example, at t — 120 nm and an 
angle of incidence (rotation angle) of 6' = 21° the ratio d±/s ± 
is twice as large as d/ sq (normal incidence) while the diffrac- 
tion angles, to good approximation, only increase by the ratio 
d/d±. This effect increases the experimental resolution. Due 
to the shadowing, however, the fundamental results obtained 
earlier for atom diffraction at normal incidence cannot be car- 
ried over unchanged. 

While we illustrate our results using the experimentally 
most interesting case of helium trimer diffraction, the gen- 
eral findings of this work equally apply to other weakly-bound 
trimers, possibly consisting of non-identical atoms. The arti- 
cle is structured as follows. In Section [H] we derive, from 
quantum mechanical scattering theory, the transition ampli- 
tude for an atom diffracted from a bar of finite thickness. In 
Section|n]we construct a periodic transmission grating from 
many bars and introduce the notion of a slit at non-normal 
incidence. We show that if the slits are not aligned with 
the direction of periodicity a characteristic asymmetry of the 
diffraction pattern arises which went unnoticed in a previous 
experiment |35]. The asymmetry is relevant for the precise 
evaluation of experimental data. In Section HVl we outline the 
general scattering theory approach to trimer diffraction and 
work it out in Section [V] for non-normal trimer diffraction 
from a grating. In Section IvTl we provide the link between 
diffraction data and the trimer bond length and discuss aspects 
of a helium trimer diffraction experiment. 

II. ATOM DIFFRACTION FROM A DEEP BAR 

In a typical beam diffraction experiment with an average 
beam velocity of the order of v = 500 m/s the kinetic energy 



per atom is a few tens of meV, much less than electronic ex- 
citation energies of the atom. Therefore, we treat atoms as 
point particles and neglect their electronic degrees of freedom 
throughout this article. The de Broglie wavelength A^b as- 
sociated with the atomic motion is typically of the order of 
0.1 nm whereas a typical length scale of the scattering object 
is d — 100 nm. We shall, in the following, refer to this relation 
as the diffraction condition: 

A dB <g:d . (1) 

The free Hamilton operator for an atom of mass mis Ho = 
p 2 /2ra. The interaction between the diffracting object and the 
atom will be described by a Lennard-Jones |36] type surface 
potential W(x) where x is the position of the atom. This inter- 
action exhibits a strongly repulsive core at a distance I from 
the diffracting object of the order of the atomic diameter, and 
it passes into a weak attractive -C3// 3 van der Waals potential 
at I > 1 nm 1121 13711 . Due to the low kinetic energy of the 
atoms in the beam it is sufficient for the purposes of this work 
to model the repulsive part of the interaction by a hard core. 
The attractive part will be omitted for the moment and will be 
be included later in Sec. 11111 

Generally, the scattering state |p', +) for an atom with inci- 
dent momentum p' and positive energy E' = |p'| 2 /2m satisfies 
the Lippmann-Schwinger equation |38] 

|p\+> = |p'> + G (£' + ;0)W|p',+> (2) 

where Go(z) = [z - #o] 1 is the free Green's function, or re- 
solvent. Denoting the atom transition amplitude associated 
with the potential W by 

? at (p;p') = <p|W|p',+>, (3) 

where p = p' + Ap is the outgoing momentum and E = 
2m, the S matrix element has the usual decomposition 

<p|S |p') = 5 (3, (p - p') - 2ni6(E - E') f at (p; p'). (4) 

In many applications the diffraction object may be effec- 
tively regarded as translationally invariant along one direction 
whence the diffraction process can be treated in two space di- 
mensions. This is the case, for instance, for diffraction from a 
slit if the vertical spread of the focussed atom beam is much 
less than the physical height of the slit. In this article we 
shall always assume the scattering object to be translation- 
ally invariant along the X3-axis. To adapt the notation, we de- 
note the two-dimensional projections into the (xi,x<i) plane 
of all three-dimensional vectors, such as p, by their corre- 
sponding italic letters, such as p. As scattering does not occur 
along the X3-axis a delta-function £(A/?3) expressing momen- 
tum conservation can be extracted from Eq. Q, leaving a two- 
dimensional transition amplitude ? at(2, (p;p') which satisfies 

f t (p;p') = 6(Ap 3 )f t( - 2) (jp;p'). (5) 

In order to derive an expression for f at(2) (p;p') we project 
Eq. (|2ji into configuration space. The full wave function 
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tfr(k',x) = 2nh{x\p' , +), where k' = p'/ft, is a sum of the 
incident part i^m C (k', x) - 2nti{x\p') and the scattered part 

<Ascatt(fc', x) = 2nh(x\Gf } {E' & + id) W\p', +). (6) 

Here, £' (2) = \p'\ 2 /2m, and the two-dimensional Green's func- 
tion, or resolvent, is Gq 2) (z) = [z - p 2 /2m]~ 1 . 

If the scale of an object is large compared to the wavelength 
of visible light it is well known that the diffraction about the 
forward direction depends only on its (two-dimensional) sil- 
houette as seen from the direction of the illuminating light, 
e.g. as for a disk and a ball. In two dimensions the silhou- 
ette of a diffraction object is simply a straight line, here called 
shadow line (line J{ in Fig. [0 and, as in optics, due to the 
diffraction condition Q the scattered part i(f scatt of the wave 
function can be approximated at small scattering angles about 
the forward direction [39]. Neglecting the attractive part of 
the potential W(x) the repulsive hard core imposes Dirichlet 
boundary conditions on the circumference of the diffracting 
object. Denoting the Green's function in configuration space 
by 



shadow line, and p' is the parallel component (cf. Fig.Q. In 
accordance with the Babinet principle, the transition ampli- 
tude shows the characteristic behavior of an optical slit func- 
tion. The Babinet principle of wave optics states that two 
complementary objects, such as a slit and a bar of the same 
width, cause the same diffraction pattern outside the direction 
of illumination (forward scattering) lf33tl . 



P 




Gf(\k'\-x,x') 



*L(x\Gf(E'W + iQ)\x>) (7) 



and using the Green theorem one finds, after some algebra, 



<Ascatt(fc 



\x) =i I dai 



Gf{\k'\-x,a)^-^ mc (k',a) 



d 



ifr- mc (k', a)—G 'Qk'\; x, a) 
da] 



(8) 



Here, da2 is the infinitesimal line element along M, and d/dai 
denotes the normal derivative (cf. Fig.Q. The Green's func 
tion can be expressed in terms of a Hankel function ET 
chap. 3.10]: 



G%>Qk'\;x,x') = -H£ > (\k'\\x-x'\). 
Using the asymptotic expansion |4l] chap. 9.2] 

H^{\k'\\x-x'\) W ~ c 



(9) 



-i\h'\x'-xj\x\ 



,i(\k'\\ X \-7T/4) 



n\k'\ 



VN 



(10) 

the far field (Fraunhofer) limit of tAscatt(fc', x) can readily be 
calculated. Inserting the expansion into Eq. ^ shows that the 
vector k = \k'\x/\x\ in the first exponential should be iden- 
tified with the outgoing wave vector, and p = tik with the 
outgoing momentum. Comparing this expression with the far 
field limit of Eq. (J6J one arrives at the two-dimensional tran- 
sition amplitude 

f at(2 W0 = 4S^4r C da ie -^l* (Ha) 
2 {2nYmh J_ A/2 



I -A/2 

i + Pa, sm(Ap a2 A/2H) 
"2 (2n) 2 mh Ap a J2h ' 



(lib) 



where A denotes the length of the shadow line M. Further- 
more, p' is the momentum component of p' normal to the 



FIG. 1: Two-dimensional geometry of atom diffraction from a 
diffraction object (bar) of finite thickness. The incident momentum 
of the atom is denoted by p'. The straight shadow line M, which 
plays the role of the silhouette, divides the bar into an "illuminated" 
part and a "shadowed" part. Also drawn is the adapted coordinate 
system («i , a^) with the «2-axis centered along the shadow line. 



III. ATOM DIFFRACTION FROM A DEEP GRATING 
A. The atom slit function of the deep grating 

In diffraction experiments one often employs transmission 
gratings to enhance the measurable diffraction peak intensities 
by a factor of A^ 2 , where is the number of coherently illu- 
minated bars. One grating bar is simply a special case of the 
general scattering object considered in the previous section. 
In the following we arrange identical bars to create a reg- 
ularly spaced periodic transmission grating. While the sim- 
plest and most familiar situation in which the bars are aligned 
along their common shadow line («2-axis, cf. Fig.[TJ is treated 
in virtually every textbook on optics (e.g. Ref. IB3I0 . we are 
not aware of a more general treatment where the individual 
shadow lines are not parallel to the alignment axis. This situ- 
ation arises naturally, however, for diffraction from a grating 
with bars of finite thickness at non-zero angle of incidence 
(cf. Fig. . We note that apart from the period an additional 
length scale of the grating is given by the distance between 
the bars. Depending on the angle of incidence the (projected) 
distance can become small or even zero. The diffraction con- 
dition Q, which must hold for all length scales, therefore also 
imposes a limit on the maximal angle of incidence. 

Under the diffraction condition the transition amplitude of a 
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grating of N bars with period d along the X2-axis can be writ- 
ten as the coherent sum of the spatially translated amplitudes 
of each bar: 



fglf (p; P) = J] e- i(n -^ N -^ d/n ? at(2) (p; p') (12a) 

n=0 

= H N (Ap 2 )f t(2 \p;p'). (12b) 

In the second line, the sum has been carried out and replaced 
by the grating function l !33ll 



S2 



H N (A P2 ) = 



sm(Ap 2 dN/2h) 
sin{Ap 2 d /2h) 



(13) 



whose argument is the momentum transfer along the direc- 
tion of periodicity x 2 . Eq. dl2l yields, in principle, a satisfac- 
tory description of the diffraction problem in terms of atom 
scattering from a bar. The literature on optics, however, com- 
monly adapts a complementary point of view by focusing on 
the apertures (slits) between the bars rather than on the aper- 
ture stops (bars) themselves. This is expressed, for example, 
by the Kirchhoff integral of optics [33]. In recent work this 
viewpoint has also proven very useful in the field of atom and 
molecule diffraction from a transmission grating: small quan- 
tum mechanical effects such as the van der Waals interaction 
1121 Il7ll between the atoms in the beam and the grating as 
well as the finite size of the helium dimer 1 3 ] manifest them- 
selves as an apparent reduction of the slit width of the grating. 
In these articles the respective quantities could be determined 
quite precisely by comparison of the reduced slit width with 
the true geometrical slit width. Unlike for grating diffraction 
at normal incidence, however, it is initially not evident how 
to define a slit in the case of non-normal incidence: the cor- 
rect choice may depend on the angle of incidence. We pro- 
vide, therefore, a mathematical prescription which converts 
Eq. Jl 2i into an expression in the style of the Kirchhoff inte- 
gral. This is achieved by introducing, between every pair of 
adjacent bars, a new coordinate system [s\, s%) as depicted in 
Fig. |3 such that the S2-axis meets the boundaries of the two 
bars at their respective shadow lines. The length of the re- 
sulting straight "slit line" S will be denoted by So- We now 
substitute the integration variable a 2 of Eq. (II 1 al by 



s 2 = 



So 
2 



as a 2 § 



where Ap S2 denotes the momentum transfer parallel to the s 2 - 
axis. Similarly, we denote the components of the incident 
and outgoing momenta with respect to the si-axis by p' s and 
p si , respectively. (An explicit expression for these momen- 
tum components in terms of the geometry of the grating will 
be given below in Eqs. (I27> and (12 8 1 for the transmission grat- 
ing of Fig.|4]) Using the identities 



Ap a2 A + Ap S2 So 

Pa, + Pa: + Ps, 



Ap a 



Ap S2 



Ap 2 d , 

Pi + p'l 

Ap 2 



(14a) 
(14b) 



which hold because of energy conservation, and the abbrevi- 
ation D = dAp 2 / Ap S2 , the transition amplitude of the single 



So/2 



\X2t 






S^Aso/2 

/ . 











3K 



d 



FIG. 2: Geometry of atom diffraction from a grating of identical 
deep bars with period d along the xi direction. The slit line <S, which 
connects the shadow lines !R of adjacent bars, provides a general- 
ization of the notion of a slit at normal incidence. It lends itself to 
the application of the Huygens principle. Also drawn is the adapted 
coordinate system {s\ , s 2 ) with the ,s 2 -axis centered along the slit line 
and the angle a by which the coordinate systems (x t ,X2) and (s\, S2) 
are rotated with respect to each other. 



bar Eq. dl 1 al l becomes 



2 (2n) 2 mh 

S I2 

ds 2 e - iAP °2 S2/n 

D/2 



JSo 



D/2 

12 



(15) 



Keeping in mind that the integration variable 52 was substi- 
tuted for a 2 for a single bar the following geometrical interpre- 
tation is possible: in the integral from So/2 to D/2 in Eq. (I15l l 
the variable s 2 represents the position along the upper half of 
the slit line S on one side of the bar; similarly, in the inte- 
gral running from -D/2 to -So/2, s 2 is the position along the 
lower half of the next slit line at the other side the bar. Insert- 
ing Eq. ( 1151 into Eq. dl2al i the half slit lines of the adjacent 
bars can be joined to yield N — 1 slits between them. Collect- 
ing all terms and introducing a "slit function", 



r D/2 

a at (p';Ap S2 )= ds 2 ex P (-iAp S2 s 2 /n) T at (p';s 2 ), (16) 

J-D/2 



where the transmission function T at (p';s2) inserted here is 
unity inside the slit >S and zero otherwise, the transition am- 
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plitude of the grating can be written, for AT » 1, as 



4£W)- iPsi+P ' s > 



2 (2n) 2 mfi 
sin (Ap 2 Nd/2h) 



(17) 



H N .i(Ap 2 ) a at (p';Ap S2 ) 



Ap s J2h 



The first term in curly brackets is sharply peaked about the 
forward direction, and in the limit N — > oo, using Eq. dl4bt . it 
simply reduces to 2nhy-6(Ap2). The second term, which is a 
product of the grating function Jl 3i and the slit function dl 6I >. 
generates the familiar diffraction pattern of a grating ll34ll . The 
n-th order principle diffraction maximum appears at the mo- 
mentum transfer 

n2nh 

Introducing the angle of incidence ff such that p' { = \p'\ cos ff 
and p' 2 - \p'\ sin ff (cf. Fig.|2j, and equivalently the diffraction 
angle 9 such that p\ - |p|cos6i and p 2 - \p\sind, the n-th 
order is located at the angle 8 = 0„ satisfying 



sin 6„ = sin ff + 



n2nh 



(18) 



Generally, the diffraction intensities are proportional to the 
scattering matrix element |(p|5 |p')l 2 where the components of 
the outgoing momentum p must be evaluated at the angle 6„. 
Inserting Eq. Mil into Eq. @ one finds, after some algebra, 



P^+PmX \a dt (p';Ap S2 )\ 
2p> J kr*(p';0)| 2 



(19) 



where the intensity 7q of the zeroth diffraction order serves as a 
normalization constant depending on the experimental count- 
ing rate and where, from Eq. ( I16> . |a at (p'\ 0) | 2 = Si (in the 
absence of the van der Waals interaction considered below). 

We now discuss the intensity formula Eq. ( I19> to explain 
the origin of the asymmetry of the diffraction pattern. The 
slit function a' dt (p'\ Ap S2 ) is an even function of the mo- 
mentum transfer component Ap Sl . The geometrical factor 
(Psi + / , j,) 2 /(2p' Jl So) 2 (the component p Sj depends on Ap Sj 
through conservation of energy and momentum) can be shown 
to introduce, for positive incident angle 6', a slight attenuation 
of the slit function at positive Ap S2 and likewise an intensifi- 
cation at negative Ap Sl . The product of the slit function and 
the geometrical factor serves in Eq. ( I19> as an envelope func- 
tion which is probed by the grating function at the momen- 
tum transfer Ap 2 rather than at Ap Sl . Since Ap Sl and Ap 2 
are not proportional to each other this probing is not symmet- 
ric for positive and negative diffraction angles. This is de- 
picted in Fig. |3 and it leads to the characteristic asymmetry 
of the diffraction pattern of a deep grating at non-normal inci- 
dence. Expanding Ap s „ using Eq. ( I14bi . into a power series 
in (Ap2/ cos ff) through second order, 



cos(ff + 9') cos 6' 



Ap 2 + tan6»' - tan(o- + ff) I Ap 2 



2\p'\ 



\cos ff 




+0.02 +0.04 +0.06 +0.08 

Momentum transfer Ap S2 / h [ 10 nm ] 

FIG. 3: Asymmetric probing of the envelope function by the grat- 
ing function at non-normal incidence. Both curves show the product 
(p.,, + p' I1 ) 2 /(2pj 1 5o) 2 |a al (p'\ Ap„) | 2 as it appears in Eq. fHJl ver- 
sus the momentum transfer Ap S2 . The solid curve refers to the pos- 
itive values of Ap S2 on the horizontal axis. The dashed curve refers 
to the negative values of Ap„ which have been mirrored onto the 
positive axis for comparison. The circles and triangles on top of 
the solid and the dashed curves, respectively, mark those values of 
Ap S2 where the grating function H N (Ap2) probes the slit function, 
i.e. where Ap 2 = n2nh/d is satisfied for n = 0, ±1, ±2, .... Their 
intensities correspond to the measurable diffraction peaks. For the 
calculation of this figure the grating cross section of Fig.|5]was used 
with the beam parameters ff = 21° and \p'\/Ti = 10 mrr 1 . 



(20) 



the leading non-linear term is seen to vanish like Ap2/\p'\ rel- 
ative to the linear term. Accordingly, the asymmetry is less 
pronounced for smaller diffraction orders, for faster beams 
and, in the case of molecules, for heavier molecules. While 
at \p'\/ti = 10 nm -1 (corresponding to a 4 He beam at v ~ 
160 m/s), as seen in Fig. the quadratic term in Eq. J20l > 
is responsible for a ±8% deviation of the positive and neg- 
ative fifth diffraction orders, respectively, its contribution re- 
duces to +0.7% at v = 1800 m/s. This smallness explains 
why the asymmetry has previously been missed (cf. Fig. 5 in 
Ref. 1 35]). Clearly, in the thin grating limit a — > the mirror 
symmetry is recovered in Eq. OOb . 

In the derivation so far no comment has been made about 
the inclusion of the attractive van der Waals interaction be- 
tween the atom and the grating. It can be accounted for 
through the transmission function r at (p'; ^2) in the slit func- 
tion \\6\ as outlined in Appendix lAl and in Refs. H5T1 . 
Unlike the case of normal incidence, at non-normal incidence 
the influence of the van der Waals interaction may be differ- 
ent on each side of the slit, introducing, in principle, an addi- 
tional source of asymmetry to the diffraction pattern. Numer- 
ical comparisons using the explicit expression of AppendixlAl 
for the transmission function demonstrate, however, that this 
effect is minor. 
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B. The atom diffraction pattern of the deep grating 

The quantitative evaluation of experimental diffraction data 
requires to determine a set of parameters describing the geom- 
etry of the particular grating (cf. Fig.0} as well as the van der 
Waals interaction coefficient C3 (cf. appendix [A}. Previous 
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FIG. 4: Geometrical cross section of a custom diffraction grating 
as it has been used in matter diffraction experiments H [121 . The 
parameters and their typical values are: period d = 100 nm, slit width 
s = 60 nm, thickness t = 120nm, and wedge angle /? = 6°. The 
angle by which the coordinate systems (x\ , x{) and (s[ , S2) are rotated 
with respect to each other satisfies cot a = tan/? + s /t and 5 sin a = 
t. At the above parameters a « 58°. The characteristic shape of the 
bars is reminiscent of the lithographic production process II ill . 

work has shown that an immediate numerical fit of Eq. dl9> 
to experimental data does not reliably determine these param- 
eters. Analogous to the procedure developed in Ref. 11211 for 
diffraction at normal incidence we therefore introduce a two- 
term cumulant approximation of the slit function. To this end 
we rewrite Eq. d!6l > as 



fl a V;Ap S2 ) 



iAp s 



JAp !2 S /2ti q 



&p s 



_ e -i&p S2 S /2h 



+ / A P»2 



(21) 



n ) \ h 

where the functions <& + (k) and are defined by 

Here, r at '(p'; 52) denotes the derivative of the transmission 
function with respect to its position argument. As (f) ± (0) = 1 
the logarithm of ± (/c) can be expanded into a power series of 
the form 



ln<t ± ( 



Z(±iK) J , 
-]T R 1- 

7=1 



(23) 



The complex numbers /?* and Ri, which are known as the 



d23l . One finds 



and 

R -(yr-< s f> 2 - 2 r" f ^'( p ' ;i (f- f ))- 

(25) 

Using the explicit form of the transmission function derived in 
the appendix the length scale of the cumulants can be shown 
to be set by the parameter ylCsl(ftv). For helium and a SiN x 
grating C3 ~ 0.1 meV nm 3 11211 . Therefore, for the purposes 
of this work it is sufficient to truncate the expansion Eq. J23l > 
after the second order. Inserting the first two terms into the 
slit function Eq. d_ 1 1 > and introducing the four quantities 



S eB = S - Re(fl|+/?7), 



A = lm(R+ +R7), 



r=Im(fl[-flp, S= y^Re(^ +R-), 

the n-th order diffraction intensity relative to the zeroth order 
is given, within this approximation, by 



In 
k 



2 K yfii 



+ A 2 



exp [-(Ap. !2 2:) 2 //z 2 ] exp[-FAp. !2 /^] 



sin 2 [A/?. S2 5 e ff/2^] + sinh 2 [Ap j2 A/2ft] 
[Ap S2 /2ftf 



(26) 



Here, the momentum components are to be taken explicitly at 
the angles 



p Sl = cos(a + Q n ) 
p' Si = \p'\ cos(a + 6»') 

and the momentum transfer is given by 

Ap Sl = \p'\ [s'm(a + 6„) - sin(a + 6')] 



(27a) 
(27b) 



(28) 



first two cumulants, are uniquely determined by Eqs. \22\ and 



where a denotes the angle shown in Fig. 0] by which the 
coordinate system (s\,S2) is rotated with respect to {x\,X2). 
The diffraction intensity formula J26I is, though more com- 
plicated, reminiscent of the result for normal incidence de- 
rived in Ref. Ill 211 : The long fraction involving the sine of S e ff 
resembles the Kirchhoff slit function for a slit of "effective 
width" S eff whose intensity zeros, though, are removed by the 
hyperbolic sine function involving A. The Gaussian exponen- 
tial reflects the suppression of higher diffraction orders. The 
asymmetry /„ + /_„ of the diffraction pattern is now embod- 
ied by the asymmetry in Ap Sl as a function of n in Eq. d28t >. 
The second exponential in Eq. \26\ accounts for the minor 
additional asymmetry in the diffraction pattern due to the dif- 
ferent influence of the van der Waals interaction on both sides 
of the bar. A comparison of the diffraction intensities calcu- 
lated from Eq. \\% and the approximation \26\ is displayed in 
Fig.0 




-2 2 
Diffraction order n 

FIG. 5: Diffraction intensities of a Helium atom beam at v = 500 m/s 
from a d = 100 nm transmission grating as displayed in Fig. [4] at 
ff = 21° angle of incidence. The solid curve was calculated using 
Eq. 1191 . the dashed curve shows the two-term cumulant approxima- 
tion I26i . To guide the eye these functions are shown continuously. 
The circles on top of the solid curve at integer n mark the experimen- 
tally accessible diffraction orders I„/Iq. 



FIG. 6: One of three possible sets of Jacobi coordinates. The vector 
p (1) points from the center of mass of the subsystem (23) to atom 
1. The vector r (23) is the relative coordinate of the subsystem. The 
coordinate R (not shown) corresponds to the center of mass position 
and is, therefore, identical for all three sets. 
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J 





(30) 



By a numerical fit of Eq. J26i to experimental diffraction 
data the effective slit width S e e, amongst the other param- 
eters, can be determined accurately and allows for further 
comparison between theory and experiment along the lines of 
Ref. fl2il . For completeness, we note that 



S eff = So - Re f S0/ ds 2 [l - r a V; s 2 )l . 

J-Sa/2 



(29) 



This means that the geometrical slit width So appears to be re- 
duced by the average deviation from unity of the transmission 
function r at (p'; S2). 



where 1 and denote the 3x3 unit and zero matrix, respec- 
tively, and M = m\ + ni2 + mi is the total mass. It is sufficient 
to restrict the combinations of indices to the ascending permu- 
tations 



(ijk) = (123), (231), (312). 



(31) 



The transformation between different sets of Jacobi coordi- 
nates can be derived from Eq. J30i . It takes the form 



(32) 



( R > 




( R 1 


p U) 


= j ui) 








{ r« J 



IV. TRIMER DIFFRACTION THEORY 

A trimer, in the scope of this article, is a three-atomic 
molecule which is weakly bound by pair interactions 1 42 ] . 
Again, the atoms themselves are treated as point particles. 
An additional three-body interaction between the atoms is as- 
sumed to be negligible. Central to later applications will be 
the helium trimer 4 He3 in which case at least these assump- 
tions are expected to be valid l3(ill . 



A. Three-body bound states 

The masses of the three atoms at positions r,, for i = 1, 2, 3, 
will be denoted by m, and are assumed to be of the same order 
of magnitude. The interaction between atom j and k is mod- 
eled by a potential Vjtdr^'l) where r ( -*' = r^-r^ is the relative 
coordinate. We introduce the Jacobi coordinates R,p (!) ,r^ 
sketched in Fig. |6j which can be expressed in block matrix 



where the block matrix JT^ is given by 



(ji) 



i 1 






^ 
1 



rtij +m* (in . +nij c )(m^ +m,- ) 

-1 — ?-i , 



(33) 



The three matrices J'^ satisfy the relations det^ <J,) = 1 
and J"0')j-('<:)jj-(<:;) _ \ Expressed in Jacobi coordinates the 
Hamilton operator for a free trimer is given by Ho + V where 



1 -2 M 
H = P + 

2M 2mi(mj + nik) 



&2 + rrij+mk uk)2 
2mjink 



+ Vki 



2 j« 



rtij + mk 



(34) 



Here, P, q (,) , and are the conjugate momenta associated 
with R,p w , and r^, respectively. Denoting the eigenstates 
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of the center of mass momentum P by |P) the full trimer states 
can be written in product form |P, <f> 7 ) = |P)|0 y ) satisfying 



[H + V]\P,<p 7 ) = E\P,<p y ) 
with energy eigenvalues 

|P| 2 

E = —+Ey 

2M 7 



(35) 



(36) 



where E y is the negative binding energy of the trimer bound 
state \<p) Y . 

The representation of a trimer state by its wave function 
depends on the particular set of Jacobi coordinates. Denoting 
the common eigenstates of the relative momentum operators 
q w and p'-* by |q, p)i*, where q and p are the corresponding 
eigenvalues, we introduce momentum space wave functions 
by 



^* } (q.p) = U*<q>Pl0y>- 



(37) 



Because of the transformation Eq. fl32i wave functions with 
respect to different sets of Jacobi coordinates can be chosen to 
satisfy the transformation relation 



»f ! V\p ( *°). 



(38) 



Since the discovery of the Efimov effect [ 3 1 ] in 1 970 the he- 
lium trimer 4 He3 has received much attention |4, 5, 6, 20, 21, 

predicted to possess, apart from its quite tightly bound ground 
state (Eg = -126 mK), a weakly bound Efimov-type excited 
state (E e = -2.3 mK) |4]. As both states have zero total angu- 
lar momentum the corresponding wave functions only depend 
on three coordinates which may be taken asp = \p\,r = |r|, 
and the angle between p and r. A common way to visualize 
trimer wave functions is to draw the hyperspherical probabil- 
ity density P(R): The hyperradius R, which is independent 
of the choice of the set of Jacobi coordinates, is defined as 
poR(p, r) 2 = \mp 2 + jinr 2 and the corresponding probability 
density can be calculated according to 



P(R) 



-s 



d 3 pd 3 r\cf> y (p, r)\ 2 6(R(p, r) - R). 



(39) 



The purpose of the "mass" parameter po is to ensure that the 
unit of R is length. The numerical value of po is, in principle, 
arbitrary as it simply scales R. Fig.Qdisplays the hyperradial 
probability densities for the two helium trimer states, using 
po/m = |, which were calculated numerically from the mo- 
mentum space Faddeev equations 1441 14511 in the unitary pole 
approximation 1 46] based on the TTY potential |47]. 



The corresponding configuration space wave functions 
4> y ^ k \p, r) are defined analogously. In order not to overload 
the notation we will omit in the following, where possible, the 
indices of the relative coordinates: if not denoted otherwise 
we implicitly use (ijk) = (123) whence p = p <l \ r = r (23) and 
q = qW p = p< 23 >. 




He ground state 



He 3 excited state 



\x 10 



100 



Hyperradius R [nm] 



FIG. 7: Hyperradial probability densities P(R) according to Eq. 1391 
of the two theoretically predicted bound states of the helium trimer 
4 He3. The states were calculated numerically using the momen- 
tum space Faddeev equations 144. 14511 in the unitary pole approxi- 
mation [46] based on the TTY potential 1 47] for the helium-helium 
interaction. Clearly, the excited state, with its expectation value of 
the hyperradius of -Jpojm (R) c = 10.1 nm, is spatially more ex- 
tended by almost one order of magnitude than the ground state with 
^l-io/m (R)„ = 1.1 nm. The scaling of the horizontal axis is logarith- 
mic. 



B. Scattering theory approach to trimer diffraction 

We now proceed to the diffraction of a trimer from an ex- 
ternal potential 

W(r,,r 2 ,r 3 ) = WiQei) + W 2 (r 2 ) + W 3 (r 3 ) (40) 

where the W,(r,) are the interactions of the individual atoms 
with the diffraction object. The full Hamilton operator is given 
by 



H = H + V + W. 



(41) 



By virtue of this structure, which is formally identical to that 
of dimer diffraction fl6ll . we may carry over the fundamen- 
tal algebraic relations from previous work. Introducing, as in 
Ref. 1 16], the resolvents 



(42a) 
(42b) 
(42c) 



G Q (z) = 


lz- 


HoV 1 


Gv(z) = 


[z- 


Ho - vr 


G w (z) = 


[z- 


H - wy 



and the two-body T matrices in three-body space, 

T V ( Z ) = V + VG v (z)V 
T w (z) = W + WG w (z)W 



(43) 
(44) 



for the potentials V and W, an AGS type 1 52] transition oper- 
ator Uvv can be derived which satisfies the equation 



Uyy — Tw + TwGqTvGoUw- 



(45) 
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In particular, the transition amplitude is given by the matrix 
element of this transition operator 11611 . 

f(P, <p y ; P', 4> Y ) = <P, <p y \Uvv(E' + /0)|P', 4> Y ), (46) 

and it determines the S matrix associated with H as 



<P,0JS w |P'A<> 



: (5 (3) (P - Y')5 y y 

-2md(E-E')t(P,(f> y ;P',</>y). (47) 



W(i"i,r2,r3). Splitting W according to Eq. J40i into the in- 
dividual potentials W,(r,) and using the Lippmann-Schwinger 
Eq. (0 the matrix element d52l can be expressed by the known 
atom transition amplitudes (0: 

(Pi.Pa. Psl Wi . Pi- P 3 . +> = <f (Pi; Pi) <5 0) (P2 " P 2 ) 
x <5 <3 >(p 3 - p 3 ) + ff( Pl ;p;) 5 (3) (p 2 - p 2 ) l?(p 3 ;p£) 
1 1 



We shall in the following impose the condition 



J'|2 



E'-E 3 + iO 



x f(p 3 ;p 3 ) 



/o 



i 



i 



\Ey\,\{<f>y\V\cf>y)\ ^ 



2M' 2M 



(48) 



E'-E 2 + iO E'-E 3 + iO 



+ ff(Pi;p'i)f2 t (p 2 ;p 2 ) 

+ cycl. perm. (53) 



which ensures that the internal energies of the trimer (both 
binding energy and potential energy) are much smaller than 
the external energy associated with the center of mass mo- 
tion. For a helium trimer beam at an incident beam velocity 
of v = 500 m/s, for example, the center of mass kinetic en- 
ergy \P'\ 2 /(2M) w 16meV (corresponding to 180 K) exceeds 
the trimer ground state energy by more than three orders of 
magnitude. Using the Schrodinger equation for bound trimer 
wave functions y (q, p) the condition J48i can be shown to 
entail the relations 



Here, "cycl. perm." indicates that all explicitly shown terms 
on the right hand side of Eq. d53l > should be repeated with their 
indices permuted once and twice, in ascending order. Apply- 
ing again the condition of the weak binding energy J48I the 
complex energy denominators can be approximated. Firstly, 
using the principal value formula (x + i0)~ l = —mS(x) + fx 1 
it is possible to approximate 



£ (3, (P2 - P 2 ) 



E' 3 -E 3 + iO 
£ <3) (P2-P2) [-2ni6(E 



E\ -Ex + iO 



E[)] (54) 



|q|«|P|,|P'| and |p|«|P|,|P'| 



(49) 



These state that the wave functions of the trimer are concen- 
trated in momentum space at relative momenta far smaller 
than the center of mass momentum. 

Under the conditions ( I48> and ( 1491 an approximation of the 
equation for Uvv (I45> to lowest order is possible and sufficient 
II53L chap. 3.4] whence the transition amplitude becomes 



where a small correction term 0(E y /E) was neglected. Sec- 
ondly, for three variables x, y, z with x + y + z — the distribu- 
tion identity 



1 



1 



1 



1 



1 



1 



x + iOy + iO 
(2;r) 2 



y + iO z + iO z + iO x + /O 
[6(x)6(y) + 6(y)6(z) + S(z)6(x)] (55) 



f(P, 



; P', 4> Y ) - <P, 4> y \T w (E' + /0)|P', <p y ). (50) can be shown to hold. If x + y + z * 0, such as in Eq. (|53|l 



We note that within this approximation the trimer binding po- 
tential V is only implicitly contained through the bound states 
\<p) y and \<p)y. The evaluation of the right hand side of Eq. (I50i 
is nontrivial. A series of approximations, all accurate within 
the condition ( I48> . may however be applied to simplify the 
transition amplitude. As the first step, the matrix element of 
Tw(E' + iO) can be shown to vary slowly under a variation 
of E' on the scale of the binding energies Ey. This allows 
to replace the energy argument of 7V in Eq. ( 15 0> by the sum 
E' y + E' 2 + E' 3 where E\ = |pJ| 2 /2;«, are the energies of the 
free atoms. Introducing two complete sets of states Eq. (I50i 
becomes 

r(P,0 r ;P',0 r ) a j d i qd i pd i qd i p' 0*(q,p)0 r (q',p') 

x <p 1 ,p 2 ,p 3 |r w (£' 1 +^ + ^+/0)|p' I ,p 2 ,p3>. (51) 

Using Eqs. (1421 and (1441 the algebraic relation Tw(z)Gq(z) - 
WGw(z) can be shown to hold. Inserting this relation the ma- 
trix element of Tw{E[ + E' 2 + E' 3 + iO) is replaced by 



for x = E', 



FJ 

C 2 



Ez, and z 



£3, Eq. J551 is 



still applicable within the same range of validity as Eq. < !54b . 
Combining the steps one arrives at the following approximate 
expression for the matrix element of TwiE^ + E' 2 + E'„ + iO): 

<pi,p 2 ,P3I7V(£; + e' 2 + e' 3 + /o)ip;,p 2 ,p 3 ) - tf ( Pl ;p;) 



x 5 (3) (p 2 - p 2 ) 6S»(p 3 - p' 3 ) - 2ni 8 (E 1 - E[) tf ( Pl ; p\) 



r(3) 



x 8 



■(3) 



(Pa - p 2 ) f?(p 3 ;p 3 ) - ^-6 (E2 - E' 2 ) 5 (E 3 - E' 3 ) 



x ff(Pi; pi) 4'(Pi; p 2 ) ff(p 3 ; p 3 ) + c y c1 - P erm - 



(56) 



Upon insertion of Eq. d56i into the trimer transition ampli- 
tude J5 1 1 several momentum integrals can be carried out by 
virtue of the momentum delta-functions. Moreover, the en- 
ergy delta-functions allow the integration of a further momen- 
tum component each. As an example, we consider the delta- 
function c5(£i - E',) in the second term of Eq. ( I56l l. Switching 
back to Jacobi coordinates, after integrating over (5 (3) (p 2 - p 2 ) 
it becomes 



<p 1 ,p 2 ,p 3 |W|p' I ,p 2 ,p 3 ,+) 



(52) 



6(E 1 -E' 1 ) = 6 



mi 

— P + q 

M H 



where |p' p p 2 ,p 3 , +) = |p' p +)|p 2 , +)|p' v +) is the scattering 
state of three independent atoms associated with the potential 



m 2 + m 3 m 2 +m 3 
P ~ — P + q Ap 



M 



in 2 



2"\ 



(57) 
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To proceed we decompose the momentum vectors into their 
components parallel and perpendicular to the incident center 
of mass momentum P': the parallel component of q, for exam- 
ple, is denoted by q\\ and the two-dimensional perpendicular 
vector is denoted by q ± . Then, by condition d49l and using 
P^l = by definition, the delta-function (I57> can be approxi- 
mated by 



S(Ei - E[) 



2=^(^+91) 



where the term involving the factor £ m 6(A^a/d) 2 is a cor- 
rection which is small by the diffraction condition Q. Inte- 
grating, according to Eq. 15 II . the second term of Eq. ( 1561 
over dp'^, the correction involving £ is pushed into the func- 
tional arguments of the atom transition amplitudes as well as 
the trimer wave functions. In both cases, however, these func- 
tions vary slowly on the scale of %P\\ such that the correction 
can be shown to be negligible altogether. Similar simplifica- 
tions are readily derived for the other energy delta-functions 
in Eq. ( BoY 

Combining all steps the general expression for the trimer 
transition amplitude subject to the diffraction condition and 
the condition of weak binding energy is obtained. Before- 
hand, however, it is helpful to introduce two abbreviations. 
Firstly, we express the atom transition amplitude in terms of 
the momentum transfer Ap ; = p ( - pj as 

7?(p,;Ap,) = f,(p,;p;). (58) 
Secondly, we introduce a molecular "form factor" 

F yr (q;p)=J <*V<iV0;(q'.p') Mq + q'.P + P')- ( 59 ) 

With this notation the trimer transition amplitude assumes the 
form 

ni2 + mj 



r(P, <p Y ; K M « { FyY (0, -^pAP,; o)? (^P; Ap) 

^ f *P. 
m 2 P\\ J 



— 2m 



M(m 2 



■at m 3 



M 



P;AP„ 



m2 + mi 



m 2 



Fyy (0, -q ± ; 0, -p J tf (^P; AP h ^AP ± + q x ) 



x tf f^.P;0,^AP ± 



M 



n%2 

~m' 



;;?2 



X 7« plp ; ,— AP ± 

3 M " 



m.2 + 
mi + W3 



q_L + P_l | 

q ± - p ± j | + cycl. perm. 

(60) 



where again use has been made of condition ( I49> to simplify 
the functional arguments of the atom transition amplitudes 
where possible. 

V. TRIMER DIFFRACTION FROM A DEEP GRATING 



A. The trimer slit function of the deep grating 



In this section the general trimer transition amplitude will 
be evaluated for diffraction from a deep grating. Inserting, to 
this end, the expression dl7l for the atom transition amplitude 
into Eq. d60l yields a very long sum containing a total of 42 
terms which we do not spell out explicitly. These 42 terms can 
be classified by the number (zero through three) of references 
to the slit function in each. In the context of a many-body 
multiple scattering series expansion 15311 they may respec- 
tively be interpreted as forward (9 terms), single (18 terms), 
double (12 terms), and triple-scattering terms (3 terms). A 
cumbersome and little elucidating calculation involving the 
transformation properties of the form factor j59t with respect 
to different sets of Jacobi coordinates reveals that the eigh- 
teen single-scattering terms interfere almost destructively and 
that their net contribution is small by a factor 0(A^ /d) 2 com- 
pared to the forward and the triple-scattering terms. They may 
thus be neglected by the diffraction condition Q. Similarly, 
the twelve double-scattering terms contribute only to order 
0(Ads/d) and may also be neglected. 



The forward and triple-scattering terms, respectively, can 
be combined. As in the atom case in Eq. Q a delta-function 
6{APi) can be extracted from the trimer transition amplitude 
leaving 



t(P, (p y \ P', <f>y) = 6(AP 3 )t (2) (P, 4> 7 ; P', (f>y) (61) 



where bold italic letters again denote the two-dimensional 
projections of vectors into the plane perpendicular to the x 3 - 
axis. The two-dimensional trimer transition amplitude for 
diffraction from a transmission grating becomes 
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i P Sl + F It Pi 1 P% 

-2^Mn dMp-KMWrr - ^2^2 



J dq±dp ± 



Fyy(0,-q ± ,0;0,-p ± ,0) 



H N Up ]} sin & + f ^ AP X + qA cos tfjaf ( ^P'; APy sin(a + ff) + ^ Ap ± + <1±) cos ( a + °')) 

q± + Pxjcos(a + 9') 
q± - Px cos(ff + 9') 



X H N \\ —AP± 

. M m 2 +m- i 



°- + ''> fl, ki p ' : ¥ APi 



+ cycl. perm. 



1 ' M \ M m 2 + m 3 



m 2 
m 2 + 
m 3 



(62) 



Inserting now for each atom the slit functions d!6t and rewrit- 
ing the form factor d59t as a configuration space integral 



(q;p) = j 



d pd re 



-%-p+p-r)/S 



*(p, r) 4> Y ip, r) 



the integrations over dq ± and dp ± in Eq. i62\ can be carried 
out. The three grating functions give rise to a triple sum 
of which only the on-diagonal terms contribute significantly: 
they represent diffraction of all three atoms from the same 
bar; the off-diagonal terms, which correspond to diffraction 
of atoms from different bars, are negligible since the proba- 
bility for atoms to be spatially separated as far as the distance 
between two adjacent bars (lOOnm) is strongly suppressed by 
the bound state wave functions of the trimer. Collecting the 
remaining terms the trimer transition amplitude can be cast 
into the form 



t^liP,<p 7 ;P',<f> Y ) 



2 i2n) 2 Mh 



(63) 



x { 2nhj-8iAP 2 )S 7 y - H N iAP 2 )a^ Y iP'- AP Sl ) 



incorporates the complicated internal configuration of the 
trimer molecule through the bound state wave functions. In 
particular, the notation 

m 2 + in?, 

ri± = S 2 cosia + ) H — — p ± (66a) 

fill YYl"i 

r 2± = S 2 cos(a' + ff) - —pj_ + r x (66b) 

M m 2 + W3 

r 3± = S 2 cos(a + #') - — — r x (66c) 

M m 2 + m-} 

has been chosen to emphasize the geometrical meaning of the 
position arguments of the atom transmission functions: The 
quantities r, ± / cos(a + 9') can be interpreted as the positions 
of the individual atoms projected onto the slit line S while 
the integration variable S 2 represents the projected center of 
mass position R ± / cos(a + 6'). The trimer transmission func- 
tion d65i is, therefore, simply the product of the three atomic 
transmission functions averaged over the wave functions of 
the incident and the outgoing bound trimer state. This intu- 
itive result is a straightforward extension of the case of dimer 
diffraction H. 



Here, we introduced a trimer slit function by 

-D/2 



p 

a^,iP';AP S2 ) = 

<J — i 



dS 2 exp i-iAP S2 S 2 /h) t^,(P';S 2 ) 



-D/2 



(64) 



where S 2 can be interpreted geometrically as the center of 
mass position of the trimer along the slit line (cf. Fig. |5J 
and where, analog to the atom case, D = dAP 2 /AP Sl . Both 
Eqs. d63l > and (I64> exhibit the same structure as their atom 
counterparts. Only the new trimer transmission function, 
which appears in the trimer slit function d64b . and which turns 
out as 

r* (P'; S 2) = J d 3 p d\ f 7 ip, r) cf> Y ip, r) 



Xif\—P'; 



x 



M cos(a + 6') ) / \M cos(a + 0') 
M ' cos(a' + 9') 



(65) 



B. The trimer diffraction pattern of the deep grating 

Thanks to the formal coincidence of Eqs. fl63i and d64b with 
their counterparts Eqs. d!7l > and dl6> of atom diffraction the 
derivation of the n-th order relative diffraction intensity lj, y 
for the incident bound state <py and the outgoing bound state 
<p y can be carried over immediately. Therefore, we write the 
trimer diffraction intensities in the form 



fry 



2P> 



jg, jP';AP S2 )\ 
i< iP';0)\ 2 



(67) 



Contrary to the atom case, however, the trimer slit function 
depends on the spatially extended trimer bound states <p y and 
<py and, therefore, in general \a l " y , (P'; 0) | 2 < Sq. 

Eq. d67l > determines the diffraction intensities of both elas- 



tic (<f)y 



') and inelastic i<p y + <py) processes. Earlier works 



on the helium trimer 1 5 1 ] as well as on van der Waals dimers 
jsill have shown that diffraction orders corresponding to in- 
elastic processes are typically less intense by five to six orders 



12 



of magnitude than those of elastic processes. They are, there- 
fore, experimentally less relevant. In the following we focus 
on elastic processes. Analogous to the procedure in Section 
fTTTI ij/ can be approximated by a two-term cumulant expan- 
sion. The cumulants /?* . now depend on the trimer state <f> 7 . 
Moreover, because of the threefold van der Waals interaction 
an additional term Q r = i lm(R2 7 -R 2y ) should be retained in 
the expansion for sufficient numerical accuracy. Taking these 
generalizations into account the n-th order relative intensity 
becomes 



In 



P„ + P' 



IP 



Si -y eff.y 



+ A? 



exp [-(AP. S2 Z r ) 2 /ft 2 - AP s S y /h] 



sin 2 [AP !2 S e ff, r /2/z] + sinh 2 [(AP j2 A y /fc + AP 2 ,Q r /^ 2 )/2] 



[AP J2 /2/z] 2 



(68) 



where, analogous to Eqs. J27i and d28i . the momentum com- 
ponents are to be evaluated at the incident angle 6' and the 
diffraction angle 9„ as 

P S1 = \P'\ cos(a + 6 n ), P' Si = \P'\ cos(a + &), 

and the momentum transfer parallel to the 52-axis is given by 

AP S2 = \P'\ [sin(a + 9 n ) - sin(ff + 6')] . 

Analogous to the atom case the effective slit width S e ff,y is 
related to the trimer transmission function (1651 by the equation 



Seff.y = So - Re f ° dS 2 [l - AP';S 2 )\ ■ 

J-SJ2 



(69) 



Fig. [3] shows elastic diffraction intensities for a beam of 
4 He3 in its ground state calculated according to Eqs. d67l and 
d68i . The asymmetry of this diffraction pattern is not as pro- 
nounced as in the atom case (Fig.|5Jl. This is due to the three- 
fold mass of the trimer which entails a three times shorter de 
Broglie wavelength. Similarly, Fig.|5]shows diffraction inten- 
sities for a beam of 4 He3 in its excited state. 

Since inelastic diffraction processes are negligible an ex- 
perimental diffraction pattern of a 4 He3 beam will in general 
be well described by an incoherent superposition of the indi- 
vidual diffraction patterns of the two bound states weighted 
by their relative population numbers in the beam. In the fol- 
lowing section we first derive the trimer size effect for a pure 
beam containing trimers in only one state. Hereafter the treat- 
ment of a mixed beam will be considered. 



VI. HOW TO DETERMINE THE TRIMER SIZE 



A. The trimer size effect 



Since the effective slit width of the trimer j69t depends, on 
the one hand, on the trimer bound state (through t*' (P'\ S2)) 




0.01 : 



-4 -2 2 4 
Diffraction order n 

FIG. 8: Diffraction intensities of a pure beam of ground state 4 He3 
at v = 500 m/s from a d = 100 nm transmission grating at 6' = 21° 
angle of incidence for the grating geometry of Fig.|4| The solid curve 
was calculated using Eq. 1671 . the dashed curve shows the two-term 
cumulant approximation 1681 . To guide the eye these functions are 
shown continuously. The circles on top of the solid curve at integer 
n mark the experimentally accessible diffraction orders I^ g /I^ g . 
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FIG. 9: As in Fig.|S|but for a pure beam of excited state 4 He3. Due 
to the larger pair distance of the excited state the effective slit width 
is smaller, resulting in a considerably wider envelope function than 
for the ground state. 



and is, on the other hand, experimentally accessible (through 
I% 7 ), it represents a link between experiment and theory. Ear- 
lier work on atom and dimer diffraction revealed that the dif- 
ference So - Seff.r carries information about both the van der 
Waals surface interaction 1 12] and the size of weakly bound 
dimers |3]. The reduction of the slit width by the dimer size 
was found to be i (r) where (r) denotes the dimer bond length. 
The subsequent evaluation of helium dimer diffraction data 
yielded the experimental result (r) = 5.2 + 0.4 nm for 4 He2 

In the following we derive the corresponding size effect for 
a trimer. To this end we explicitly insert the trimer transmis- 
sion function d65l > into Eq. d69t . By definition the atom trans- 
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mission functions T'fip'f, sq) in Eq. (1651 are zero if their po- 
sitional arguments 5,2 = r,- x /cos(ar + ff ) lie outside the slit 
interval [So/2, So/2]. This fact may be utilized to reduce 
the integration interval for the center of mass position S2- at 
fixed relative coordinates p, r the interval may be limited, for 
rj_ > 0, to 



A 2 SO 
< .5 2 < 

cos(a + 6') 2 



and, similarly, for r ± < 0, to 

So A7 



So 



<S 2 < 
2 cos(a + 8') 2 

Here, the geometrical quantities 



cos(q' + 0') 



cos(a + 0') 



(70a) 



(70b) 



A? = ±r 



1 (m2 +M3 — mi 



M 



-P± + 



m2 + m3 



Pi. 



;;?3 



m2 + m3 



(71) 



and 



1112 + m?, — m\ 



M 



-Pi. + 



/7? 2 



m2 + TO3 



Pi + 



m 2 



m 2 + m-i 



-r± 



(72) 



have been introduced. Neglecting, for the moment, the van der 
Waals interaction, all atom transmission functions are unity 
inside the reduced domain of integration J70I . In this case the 
effective slit width depends only on A* and A^. Accordingly, 
we call it the geometrical part of the effective slit width and 
denote it by 



c geom _ c 
ert.y u 



/ „ Re f d 3 pd 3 r\<t>y(p,r)\ 2 
cos(a + 8) J 

X {[At + AJ] ©(rj + [A7 + A£] 0(-rJ} 



(73) 



where 0(r ± ) is the Heaviside step function. Both integrands in 
Eq. J73l > can be simplified using the transformation properties 
of the Jacobi coordinates j32l > and the wave functions ( I38l >. 
Combining the results, the geometrical part of the effective 
slit width becomes 



g geom 
eff,y 



(31) 



,(12) 



(74) 



2cos(a + 6»') 

where the expectation values in the numerator are defined as 



In Eq. H41 the symmetric term 



131) 



+ 



,(12) 



(75) 



(76) 



represents the expectation value of the "width" (projected di- 
ameter) of the trimer perpendicular to the incident direction. 




FIG. 10: Geometrical interpretation of the effective slit width for- 
mula. Clearly, the expression Idrj, 23 '! + Ir*, 31 '! + Ir^ 12 '!) may be in- 
terpreted as the "width" (projected diameter) of the trimer perpen- 
dicular to its incident direction. Taking the expectation values with 
the bound state wave function yields the expression 1761 . The multi- 
plication by [cos(a + 0')] 1 corresponds to an orthogonal projection 
onto the slit line <S. Hence the slit width 5 appears reduced by the 
projected width of the trimer. 



Therefore, as the factor [cos(a + ff)] -1 corresponds to an or- 
thogonal projection of the perpendicular coordinates onto the 
slit line (cf. Fig. II Oi l S^ m is simply smaller than So by the 
projected width of the trimer. 

In the presence of the van der Waals interaction an addi- 
tional term accounting for the deviation from unity of 



eff,y 

the atom transmission functions arises, 
slit width is the sum 



The entire effective 



>eff,y 



eff,y 



+ s 



eff.y - 



(77) 



As the general expression for the van der Waals part S is 
long and little informative we shall not give its general form 
explicitly. 



B. The size effect for three identical Bosons 

In the remaining paragraphs of this section we focus on 
trimers of three identical Bosons, such as the 4 He3. By con- 
sequence, we denote the masses by m = m, and the equal 
projected pair distances by <|r ± |) r . The geometric part of the 
effective slit width MA\ immediately reduces to 



c geom _ c 
o — 00 — — 



3 (\r ± \) y 



eff,y 



2 COS(Q' + 0') 



(78) 



Moreover, if the spatial extent of the bound state wave func- 
tion is small compared to the slit width, the van der Waals part 
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S V n W is to very good approximation given by 



I" 



! rSo/2 

Jo 



p'. s , K-^iP 

3 ' 2 cos(a + 0') 



x r" 



> , kJ + (p i -i|rj)®(p ± -i|rj)] 



cos(a + 0') 



dS'-, 



-So/2 



P' c , 



cos(a + 0') 



> , kxl-(px + £kj)0(-px-jh.l) 

T' 2 + cos(ff + 0') 



(79) 



Within the approximation J79I it is evident that S ^ w is indeed 
zero if the atom transmission functions are unity inside the slit, 
and if the spatial extent of the trimer wave function is small 
on the scale of the slit width. Therefore, if S v £ w were but a 
small correction to the full effective slit width (I77> it could be 
neglected, and the projected trimer pair distance (\r ± \) could 
be determined using Eq. (I78> . Experience with dimer diffrac- 
tion has shown, however, that the effect of the van der Waals 
interaction can be of the same order as the pair distance |3] 
and must be accounted for. Since Eq. J79l > depends on the 
full trimer wave function it cannot be used immediately for 
the evaluation of experimental data and an approximation is 
required. The integrand in Eq. i79\ is, however, slowly vary- 
ing on the scale of the variation of |0 y (p, r)| 2 . Therefore, the 
positional arguments of the atom transmission functions can 
approximately be replaced by their expectation values. This 
approach is in analogy to Ref. |3]. An analysis of the com- 
binations of p ± and r ± , using once more the transformation 
properties of the relative coordinates (I32> . shows that these ex- 
pectation values are expressible solely in terms of (\r ± \). For 
example, 



Px 



= (\r ± \) y 



and 



(px-|kxlj©(pi-i]r x ]^ = i <kxl> r . 



(80) 



(81) 



Inserting these one finds as the final form of the van der Waals 
part of the effective slit width 
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3 2 4 cos(a + 0') 



cos(a + 9') i 
(82) 



In order to test the validity of this approximation we carried 
out a numerical analysis of the error introduced by the replace- 
ment of Eq. M9\ by Eq. (1821 : if applied to experimental data 
the approximation entails, for the two theoretically predicted 
bound states of 4 He3, a systematic overestimation of (\r ± \) by 
7% ( 4 He3 ground state), or 3% (excited state). As seen from 
Fig.^2 tne approximation is more reliable at high velocities as 
the impact of the van der Waals interaction becomes smaller. 
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FIG. 1 1 : Effective slit widths versus the beam velocity v = |P'|/M for 
the ground state and the excited state of 4 He3. The numerical results 
using the full expression Eq. 1791 are shown as solid curves and the 
approximation Eq. 1821 as dotted curves. The angle of incidence was 
taken as 8' = 21°. The approximation becomes more reliable at high 
velocities as the impact of the van der Waals interaction decreases. 
Asymptotically, for high velocities, both pairs of curves approach 
their respective upper limits 5^™™ given by Eq. 1781 . 

Theoretical studies of the helium trimer commonly state the 
expectation value of the pair distance (r) itself, where r = |r|, 
rather than a component such as (\r ± \). To link experimental 
results to these, a relation between (r) and (\r ± \) must be estab- 
lished. Both predicted 4 He3 bound states are spherically sym- 
metric (zero total angular momentum). Moreover, the two- 
body scattering matrix corresponding to the He-He potential 
is dominated by the shallow s-wave bound state pole of 4 He2, 
and higher partial waves may to good approximation be ne- 
glected 1 21]. By the Faddeev equations (e.g. Ref. 1451 ') for 
the helium trimer bound state, it is then possible to derive the 
relation 



<|rxl> = 5<r>. 



(83) 



In summary, the effective slit width (1771 depends to good ap- 
proximation only on one trimer parameter, namely the expec- 
tation value of the pair distance (r). Consequently, (r) can, in 
principle, be determined from trimer diffraction data. 



C. Experimental considerations 

Using the results of the preceding sections the improvement 
in resolution through diffraction at non-normal incidence over 
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normal incidence may be estimated. The evaluation process of 
experimental data involves two main steps: Firstly, values for 
the effective slit width S e ff t7 must be obtained by numerical 
fits of the intensity formula d68i to experimental diffraction 
patterns. Secondly, (r) is determined from a fit of Eq. Ml\ 
to the values for S e «, y . The stronger the dependence of this 
procedure on (r) the more precisely (r) can be determined. 
As (r) changes the width of the Kirchhoff-like slit function in 
Eq. d68l > a natural measure for this dependence is provided by 
the number of diffraction orders under the central maximum 
of the slit function to either side of the forward direction. This 
number, which we denote here by n c , and which we treat as a 
continuous variable, can be approximately written as 



s± 

where d ± denotes the (projected) period perpendicular to the 
beam and s± denotes the projected slit width. At the angle 
of incidence 6' the projected period is d± = dcosO'. Fur- 
thermore, neglecting for this estimate the van der Waals part, 
we insert s±_ = 5 , ^°™cos(o; + 8') for the projected slit width 
of the deep grating where S^ m will be taken from Eq. J78I . 
The relative variation of n c with (r) can then be calculated and 
becomes, to leading order in (r) /So, 

1 6n c ^ 3 1 
n c d{r) 4S cos(a + 6'Y 

In contrast, at normal incidence the right hand side would 
be 3/(4 so)- Inserting the parameters of Fig. [3] yields 
(d« c /d {r))ln c w 2.6 x lO^nnr 1 for 6' = 21° as compared to 
1.2 X lCr 2 nm -1 for normal incidence. This roughly twofold 
gain in sensitivity is expected to halve the final error bars on 
(r). 

Finally, since the population ratio of the two predicted 4 He3 
states in the nozzle beam is generally unknown the situation 
of a mixed beam must be considered. To analyze this we have 
summed diffraction patterns as shown in Figs.|S]and|5]for dif- 
ferent population ratios. Hereafter, we have used Eq. d68l to 
determine, from the summed patterns, an average effective slit 
width Seff and from this an average bond length (r). It turned 
out that the such determined (r) varies almost linearly with 
the population ratio from the ground state value of (r) (pure 
ground state beam) to the excited state value (pure excited 
state beam). Therefore, three possible outcomes of an experi- 
ment are to be expected. A value of (r) * 1 nm would be at- 
tributed to the ground state and indicate a negligible (or zero) 
population of the excited state. Equivalently, a result of about 
8 nm would doubtlessly provide evidence for the excited state 
and its large pair distance. Thirdly, a value in between these 
two would indicate that both states are present and evidence 
for the excited state would still be available. A controlled vari- 
ation of the accessible beam parameters might then allow to 
influence the population ratio in favor of either state and to 
measure the pair distance for one state with less disturbance 
by the other. 



VII. CONCLUSIONS 

Motivated by the long-standing interest in the Efimov effect 
l3lll we have studied the diffraction of weakly bound trimers 
in a typical matter optics setup. As an earlier diffraction exper- 
iment for the spatially extended helium dimer ((r) = 5.2 nm) 
Q had indicated that the resolution provided by a custom 
So = 60 nm transmission grating, at normal incidence, may 
be insufficient to resolve the helium trimer ground state ((r) = 
0.96 nm predicted 0]) it had suggested itself to use oblique 
(non-normal) incidence at a rotated transmission grating for 
reducing the projected slit width. The partial shadowing of 
the slits caused by the finite thickness of the etched mate- 
rial grating has required, however, a revision of the theory of 
atom diffraction. In particular, the familiar mirror symme- 
try encountered in diffraction patterns from normal incidence 
is lifted for non-normal incidence. This effect was visible in 
Fig. 5 of Ref. 13511 but went unnoticed. It has been traced 
back to the non-alignment of the direction of periodicity of 
the grating with the shadow lines of its bars, or, equivalently, 
its slit lines. The weak attractive van der Waals surface inter- 
action, which introduces an additional but minor asymmetry, 
has been taken into account in a way similar to the case of 
normal incidence. 

Using atom diffraction as one building block, the multi- 
channel many-body quantum mechanical scattering theory ap- 
proach of Refs. I15l ll61 has been extended to derive the con- 
stitutive formulas of trimer diffraction. While this procedure 
structurally partly parallels that of dimer diffraction it is math- 
ematically more complex due to the additional atom. The re- 
sulting equations for the trimer diffraction pattern, however, 
have been readily interpretable and provide intuitive physi- 
cal insight into diffraction of weakly bound molecules: the 
significant measurable quantity is the quantum mechanical 
expectation value of the "width" (projected diameter) of the 
trimer perpendicular to its flight direction. For identical Bo- 
son trimers, such as the helium trimer, the width is related to 
the molecular bond length by | (r). This fact can be used, in 
principle, to determine (r) from a matter diffraction experi- 
ment. 

If a transmission grating of the type used earlier by Grisenti 
et al fl2ll is rotated by 6' = 21° the projected slits appear 
approximately half as wide as the nominal slits of the grat- 
ing. This leads to an estimated doubling of the resolution, 
sufficient to determine the ground state pair distance of 4 He3. 
Moreover, should the 4 He3 Efimov state, whose pair distance 
is predicted to be larger by almost an order of magnitude 01, 
exist and should its population in the beam be significant, 
it ought to be clearly distinguishable from the ground state 
solely by its size. 
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APPENDIX A: SURFACE INTERACTION 

Earlier work has shown that at beam velocities typically en- 
countered in matter diffraction experiments the effective re- 
duction of the slit width due to both the van der Waals surface 
interaction and the finite molecular size can be of the same 
order of magnitude |3|, ll 211 . A quantitative determination of 
the atom transmission function r at (p'; 52), which was inserted 
into Eq. d!6i . is therefore necessary. As in Ref. fl2ll we use the 
eikonal approximation to write r at (p'; 52) = exp[/^(p'; S2)] for 
52 inside the slit, and r at (p'; 52) = outside. The phase func- 
tion <p(p'; s 2 ) is given by |55] 

<f(p'; s 2 ) = -m- 1 [dt W surf (s(f)), v=^, (Al) 
J m 

where the straight path of integration s(t) must be taken to run 
parallel to the direction of incidence, and to cross the slit line 
>S at the position s 2 ■ The surface interaction W sur f(a;) at a po- 
sition x between two bars is calculated from the integration 
of an attractive —Cell 6 potential of Lennard- Jones type (the 
repulsive part has already been modeled by the boundary con- 



ditions in Sec.|nJ over the volume of the bars. Carrying out all 
four integrations for the typical wedge-shaped bars shown in 
Fig.|4]the phase function can be calculated explicitly. Using 
the abbreviations 



nC(, — cos#' _ cos#' 

C 3 = a = — d, so = —so, 

6 cos(a + 6') cos(a + 6') 



it reads 



<p(p'; S2) = 



c 3 



2hv cos 2 8' cos 2 (a + 9') 
g n 2 + (g„ - dY 2 - - (£12 + d)- 2 
tan 6' +tany6 

+ (fti ~ dT 2 - - (&2 ~ dy 2 ' 
tan 8' - tan/3 



where f n = f - s 2 , £21 = f + * 2 , £12 = ^ - s 2 
L (a ie>) S 0' and £22 = f - s 2 - s () . 
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